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LOW-FREQUENCY FLUCTUATION SPECTRA AND ASSOCIATED PARTICLE 
TRANSPORT IN THE NASA LEWIS BUMPY -TORUS PLASMA 


by Chandra M. Singh, * Walter M. Krawczonek, J. Reece Roth, JaeY, Hong,"*^ 

Young C. Kim, T and Edward J. Powers^ 

Lewis Research Center 

SUMMARY 

The NASA Lewis biimpy-torus experiment employs a modified Penning discharge to 
produce and heat a plasma in a bvimpy toroidal magnetic field. The strong radial elec- 
tric field associated with the Penning discharge and the strong toroidal magnetic field 
give rise to a diversity of E/B phenomena, such as rotating waves and spokes, which in 
turn manifest themselves as space-time fluctuations of the plasma density and potential. 
The ultimate objective of the work reported herein is to further understand the nature 
and origin of the fluctuations and their connection with fluctuation-induced transport. 

The approach is to monitor the density and potential fluctuations; to digitize the data; 
and to generate, with the aid of a computer, various spectral properties by means of 
the fast Fourier transform. Of particular interest is the computer -generated transport 
spectrum that indicates in a quantitative way which fluctuation spectral components con- 
tribute to transport and which do not. All experimental measurements of the spectral 
characteristics of the plasma are given in absolute units rather than relative values. 
Preliminary measurements of the transport spectrum of the ion population are given, 
and it is shown that the fluctuation- induced transport is in order-of -magnitude agreement 
with that inferred from the steady-state current flowing to the electrodes that generate 
the plasma. 


♦National Research Council - National Aeronautics and Space Administration Re- 
search Associate. 

T University of Texas, Austin, Texas. 

J Professor of Electrical Engineering, University of Texas, Austin, Texas. 
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INTRODUCTION 


The NASA Lewis bumpy-torus facility is a steady-state plasma confinement device 
in which a toroidal plasma is acted upon by stroi^ electric fields externally imposed 
along the minor radius of the plasma. It employs a modified Peimii^ discharge to pro- 
duce and heat a plasma in a bumpy toroidal magnetic field geometry produced by 
12 superconducting coils, each havii^ an inside diameter of 19 centimeters and being 
capable of producing 3 teslas on axis. The coils form a toroidal array 1. 52 meters in 
major diameter. Figure 1 is an isometric cutaway drawii^ of the facility, details of 
which are given in reference 1. A Penning discharge occurs when circular electrode 
rings are placed at the midplane of one or more of the 12 toroidal sectors. The electrode 
rings can be biased with positive or negative potentials in excess of 10 kilovolts with re- 
spect to the grounded magnet dewars. All data reported herein were taken with positive 

electrode polarity. This combination of strong, crossed electric and magnetic fields can 
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produce a steady-state deuterium plasma of density 10 to 10 particles per cubic centi- 
meter, ion kinetic temperatures of 340 to 2500 electron volts, and electron kinetic tem- 
peratures of 8 to 140 electron volts. To avoid damage to the probes, we had to restrict 



Figure 1. - Isometric cutaway drawing of NASA Lewis bumpy torus. 
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the operating conditions to the low ends of these ranges. The plasma in this investiga- 
tion was subject to a steady- state radial electric field that ranged between 50 and 200 
volts per centimeter. The results of investigations into other aspects of this plasma 
are reported in references 2 to 7. 

The biunpy-torus plasma operates in two modes (refs. 2 and 3): the low-pressure 
mode (LPM), and the high-pressure mode (HPM). The low-pressure mode is charac- 
terized, in deuterium, by an electron temperature T greater than 35 electron volts 

11 ^ 

and a neutral bacl^round density Uq below 9x10 particles per cubic centimeter 
(2. 7x10”^ torr). The high-pressure mode has T^ less than 35 electron volts and ng 
that usually are above 9x10^^ particles per cubic centimeter (2.7x10"^ torr). 

When a positive voltage is applied to the midplane electrodes, the plasma floats at 
a positive potential approachii^ that of the midplane electrode rings. The presence of 
a high positive bias on the toroidal plasma leads to a radially outward electric field 
between the plasma and the grounded magnet dewars and vacuum chamber walls. A 
radially inward electric field exists between the anode rings and the plasma. The di- 
rections of the steady-state electric and magnetic fields acting on the plasma are shown 
in a schematic of the biunpy-torus plasma in figure 2. The plasma can thus support two 
oppositely directed E/B drift motions. The boundary of the plasma is defined by the in- 
side diameter of the circular electrode rings. The plasma follows the magnetic field 

Positive 

anode 

voltage, 

r Inward-pointing electric field 
/ near anode rings 



Rgure 2. - Schematic drawing of bumpy-torus plasma with positive electrode polarity. 
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lines and narrows down as it passes through the throats of the magnetic field coils. Fig- 
ure 3 shows the plasma in the equatorial plane of the torus. 

The reservoir of free energy stored in the E/B drift mption of the toroidal plasma 
is a potential source of instabilities, fluctuations, and turbulence. Such disturbances 
can also be driven by the stroi^ radial electric fields. The physics of crossed-field 
devices is reviewed in references 8 to 10, The fluctuations in density and potential have 
been implicated in the thermalization of the directed energy of the plasma particles and 
also in the transport of particles to the electrodes (refs, 2 and 11 to 15). 



Figure 3. - Bumpy-torus plasma, viewed in equatorial plane of torus. 


Plasma transport in fusion- oriented, crossed-field devices like the NASA Lewis 
bumpy torus is complicated by the presence of large steady- state electric fields, up to 
several kilovolts per centimeter (ref. 3). Another factor that makes the study of plasma 
transport interestii^ in this device is that T. » T . (Symbols are defined in the ap- 
pendix. ) When Tj » T^, space-charge sheath formation near the electrodes is not 
expected to be governed by the Bohm criterion (ref. 16). In addition, the current flow 
to the electrodes surrounding the plasma is not ambipolar (ref. 4): Only electrons flow 
to the anodes and only ions to the cathodes. This factor, along with the large, radial 
electric fields and the nonneutral nature of the sheaths, may affect the stability and 
transport of the plasma. A study of the density and the potential fluctuation spectra may 
reveal which portions of the spectrum cause transport and which are associated with 
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thermalization of the directed energy of the plasma. Knowing the dependence (refs. 17 
and 18) of the current collected by the electrodes on plasma parameters and compariiig 
this dependence with the fluctuation-induced transport should indicate whether or not 
fluctuations are the dominant transport mechanism. 

The study of fluctuation-induced transport follows a digital approach: Density and 
potential fluctuations are monitored with the aid of probes, and the correspondii^ signals 
are digitized and stored in a computer where they are analyzed by fast-Fourier- 
transform (FFT) spectral analysis techniques. In particular, the amplitude, frequency, 
and wave numbers of the waves are determined from the computer-generated autopower 
and cross-power spectra. The FFT techniques are also used to generate the transport 
spectrum that quantitatively indicates which portions of the spectriun are contributing 
to the transport and which are not. This report deals with the low-frequency w « co , 
fluctuation spectrum and the associated fluctuation-induced radial particle transport in 
the NASA Lewis bumpy torus. The high-frequency emission spectrum of this plasma 
has been reported in reference 6. Different combinations of circular electrode rii^s of 
18 -centimeter inside diameter were placed at the midplane of the toroidal sectors. The 
numbering of the electrodes and the toroidal sectors is shown in figure 2. The fluctua- 
tion spectra and transport have been examined for a limited range of the neutral back- 
ground density and voltage applied to the electrodes. It is important to realize that this 
diagnostic method will not necessarily measure diffusive transport resulting from binary 
collisions. It does, however, measure the transport resulting from fluctuating electric 
fields, which dominate in this plasma. 


FLUCTUATION DIAGNOSTICS AND DATA HANDLING 

The day-to-day operation of the NASA Lewis bvunpy-torus experiment requires the 
acquisition of analog and digital signals of various amplitudes and time durations. The 
neutral backgroimd density, magnetic field strength, electrode voltage and current, and 
background neutral gas impurities are routinely monitored (the last by mass spectrom- 
eter). When required by the experimental program, measurements can be made of 
X-ray and neutron production, ion kinetic temperature (by a neutral particle analyzer), 
time- and space -averaged electron density (by a polarization diplexing microwave inter- 
ferometer), and 10-megahertz to 1. 2 -gigahertz electromagnetic emission (by a minia- 
ture, rf, coaxial probe). 


Data Acquisition 

This investigation uses a set of capacitive and Langmuir probes to measure the po- 
tential and density fluctuations, respectively, in the frequency range from 1 kilohertz 
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to 5 megahertz. Figure 4 is a block diagram of the data-acquisition system. The probes 
were hydraulically actuated into the plasma for short times to take data. The upper limit 
on probe dwell time was the point at which the probe tips became red hot. The maximum 
dwell time with this actuation system was 0. 7 second. Inserting the probe assembly into 
the plasma caused a change of less than 5 percent in the current drawn by the electrodes. 
We therefore assume that the gross plasma characteristics do not change significantly 
during the time the probes are in the plasma. 

The travel of the probe assembly was so adjusted that the capacitive probe in chan- 
nel 1 was located 3 millimeters inside the boundary defined by the inner circumference 
of the electrodes in adjacent sectors of the torus, or 7 centimeters from the plasma 
axis in the equatorial plane. The capacitive probe in channel 2 was 30° above the equa- 
torial plane, in the same vertical plane and at the same plasma radius as the Langmuir 
probe in channel 3. The Langmuir probe in channel 3 was in the equatorial plane at the 
same plasma radius of 7 centimeters. However, the capacitive probe in channel 1 and 
the Lai^muir probe were displaced axially 1 centimeter from each other. Figure 5 is a 
closeup view of the probes in the probe mount. With this probe arrangement (i. e. , 
channels 1 and 3) the potential and density fluctuations can be measured at the same 


Capacitive probe 
•channel 1) 



Figure 5. •’ Closeup view of probes in probe mount. 


minor radius in the plasma, without reference to the axial dimension for waves whose 
axial wavelength X|| is much larger than 1 centimeter. With the capacitive probes in 
channels 1 and 2 the azimuthal wave number k^ of the potential fluctuations can be 
measured. The probe assembly is near the midplane of sector 4, as shown in figure 2. 
The microwave interferometer is in sector 10. The axial and radial wave numbers k|| 
and k^ have not been measured in this study. Figure 6 shows the probe assembly and 
its hydraulic actuating system, and figure 7 shows this assembly mounted on the vacuum 
tank. 
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Data Processing and Analysis 


The data-processing system consists of waveform recorders, a magnetic tape re- 
corder, a microcomputer, and a computer terminal (fig. 8). In a typical experimental 
run, the probes were actuated into the plasma; the data were recorded on the waveform 
recorders; the probes were brought out of the plasma; and finally, after the sampled 
waveforms were inspected for signal saturation or other defects in the data, the data 
were recorded on magnetic tape. The tape was then sent to the University of Texas at 
Austin for processing. 

Fast-Fourier-transform (FFT) spectral analysis techniques were used to generate 
autopower and cross -amplitude spectra from which the amplitude, frequency, and wave 
number of each wave in the fluctuation spectrum were determined. The transport 
spectrum was also computed to provide quantitative insight into how much transport is 
associated with the various spectral bands. The raw waveform data were also used to 
generate such amplitude statistics as the probability density fimctions, the mean, the 
standard deviation, the skewness, the kurtosis, and higher moments of the distribution 



Figure 8. - Data-acquisition system. 
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Figure 9. - Schematic of data-processing system. 


function. The various data-processing functions performed on the taped data at the 
University of Texas are shown on the block diagram in figure 9. 


System Calibration 

Before the plasma data were processed and analyzed, the data-acquisition system 
was checked by known test waveforms from a signal generator. Waveforms of known 
amplitudes, frequency, and waveshape were applied to the input of all three waveform 
recorders. We concluded, after generating power spectra for different waveforms, 
frequencies, and sampling intervals, that there was no measurable instrument- related 
phase shift (time delay) in the data from the three channels. Further details on the 
calibration procedure and specific examples of test waveform spectra are given in ref- 
erence 5. The constant phase shift between channels discussed in reference 5 has, 
however, been eliminated in the present series of runs. 

Capacitive probes . - Capacitive probes were used to measure the low-frequency 
(w « ojp.) electrostatic potential fluctuations. It is shown in reference 19 that a ca- 
pacitive probe has a better frequency response than a high-impedance, floating Lang- 
muir probe. The attenuation and frequency response of each capacitive probe were 
measured: The frequency response was flat from 1 kilohertz to 5 megahertz. From 
the probe’s attenuation the absolute peak-to-peak amplitude of the potential fluctuations 
in the plasma could be determined. 

Langmuir probe . - The simplicity and ease of operation of a Langmuir probe make 
it a widely used diagnostic tool for the measurement of spatially resolved plasma den- 
sity. However, the theory of probe operation (ref. 20) in hot (T^ » T^), drifting plas- 
mas with strong magnetic fields is not fully understood. Therefore, a microwave inter - 
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ferometer was used to obtain the absolute value of time- and space-averaged electron 
density. A single Langmuir probe was used to measure the steady- state floating poten- 
tial and the density fluctuations. The Lai^muir probe was connected to a 100-megohm 
input impedance probe to measure the dc floatii^ potential. 

From the instantaneous floating potential the amount of dc bias from a power supply 
that is required to measure the ion or electron saturation current can be determined. 

The electron temperature was in all cases much less than the energy corresponding to 
the bias potential, hi all the runs reported herein, fluctuations in the ion saturation 
current were measured (ref. 16). A measurement of the electron saturation current 
caused the probe to draw large currents, making the probe tips glow red hot and the 
probe eventually destroy itself. All the data presented herein were taken with short 
enough dwell times to avoid any effects on the data (e. g. , noise) due to incandescence 
of the probe or arcing to the probe tip. 

The radial profile and waveform of the ion saturation current also yield information 
on the density profile and the density concentrations that reveal the presence or absence 
of rotating spokes. Electron spokes in high-vacuum magnetrons and other axisymmetric. 



(a) Background neutral gas density, oq, l.SXlO^^ particles per cubic centimeter (5.3X10-5 torr). 



2 /Asec/cm 


(b) Background neutral gas density, ng, 7X10^^ particles per cubic centimeter (2.1X10'5 torr). 


Figure 10. - Fluctuations in ion saturation current. Electrode 10 in place; anode voltage, 2 kilovolts. 
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crossed-field devices propagate azimuthally with electron density fluctuations 
n /n « 1. The spoke usually occupies only a fraction of the minor circumference. In 
high-pressure, crossed-field discharges (ref. 9), the plasma forms one or more ro- 
tating spokes, but eventually ionization, diffusion, or shock formation can lead to a 
spread in the spoke azimuth until a continuous plasma disk results. Such cases are 
characterized by nig/iig < 1.0. 

Langmuir probe measurements of the ion saturation current in the NASA Lewis 
bumpy torus show that n /n < 1. 0. A representative case is shown in figure 10. 

The lower trace is the zero-current baseline; the upper trace shows the dc, coupled 
voltage drop across a 50-ohm resistor in the Langmuir probe circuit. The shape of 
the waveform is indicative of the shape of the density distribution inside the density 
perturbation. 


EXPERIMENTAL OPERATING CONDITIONS 

This investigation was exploratory, and no attempt was made to obtain parametric 
variations of the plasma dispersion or transport characteristics. For this series of 
runs, the magnetic field was held constant at 2. 4 tesla at the coil throats and 0. 67 tesla 
at the Langmuir probe. The electrode polarity was positive, which imparts a positive 
potential to the entire plasma. This resulted in an outward-pointing electric field at 
the probe, since the probes were located at the plasma boundary, as shown in figure 2. 
Deuterium gas was used in all runs. 

Approximately 30 individual experimental runs yielded valid data. The operatir^ 
conditions of nine runs of particular interest are given in table I, including the number 
of electrode rings that were used to generate the plasma and the sectors in which the 
electrodes were located. The sector numbering convention is illustrated in figure 2. 

The probe assembly was in sector 4, diametrically opposite the sector containing the 
microwave interferometer. The positive voltages applied to the electrodes (column 5) 
can reach several tens of kilovolts, but had to be kept low during these investigations to 
avoid arcing to the Langmuir probe. 

The electrode current in table I (colmnn 6) is the total ion current flowing from the 
plasma to the power supply, and the background neutral deuterium gas density (column 7) 
was corrected for the appropriate gage factor. The electron number density (column 8) 
was measured with the polarization diplexing microwave interferometer (ref. 4), and 
the relative density fluctuation (column 9) is the ratio of one-half the peak-to-peak am- 
plitude to the mean value, as measured on oscilloscope traces. The density fluctuations 
were measured by connecting the Langmuir probe to ground through a 50-ohm resistor. 
Since the plasma was floating at high positive potentials during these rvins, the ion sat- 
uration current could be obtained without additional bias voltage. 
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TABLE I. - OPERATING CONDITIONS 


Run 

Figure 

Number of 
electrodes 

Sectors 

containing 

electrodes 

Anode 

volt- 

age, 

kV 

Anode 

cur- 

rent, 

la* 

A 

Background 
neutral gas 
density, 

particles 

cm^ 

Average 

electron 

number 

density, 

electrons 

m 

Relative 

density 

fluctua- 

tion, 

Mode of 
opera- 
tlon^ 

1 

21, 22 

11 

Not 4 

2 

0, 12 

7.0x10^^ 

1. 7x10^® 

0. 50 

LPM 

2 

11, 20 

11 

Not 4 

1 

.026 

7.0 

.46 

.20 

LPM 

3 

12, 17, 25 

2 

1,7 

1 

. 12 

7.0 

2.5 

.52 

LPM 

4 

18 

11 

Not 4 

1 

.024 

18.0 

. 50 

.40 

HPM 

5 

15 

2 

5,6 

2 

. 19 

3.6 

2.0 

00 

LPM 

6 

27, 28 

1 

10 

1 

. 12 

18.0 

6.2 

o 

CO 

HPM 

7 

16 

1 

10 

1 

. 11 

7.0 

4. 1 

1.0 

LPM 

8 

14, 19, 26 

2 

1,7 

1. 2 

.21 , 

16.0 

5. 3 

. 19 

HPM 

9 

13, 24 

2 

3,9 

1 

.082 ; 

7.0 1 

2.9 

.67 

LPM 


^LPM denotes low-pressure mode; HPM denotes high-pressure mode. 


CHARACTERISTICS OF AMPLITUDE DKTRIBUTION FUNCTION 
Theoretical Considerations 

The first step in analyzing the plasma fluctuation data was to determine the proper- 
ties of their amplitude distribution functions. In digitizing the analog signals, the wave- 
form recorders sampled 2048 equally spaced points along the waveform and at each point 
could distinguish an 8 -bit amplitude range, or 256 increments of signal amplitude. 

These 256 increments were assigned to one of 60 amplitude cells. The computer pro- 
gram places each of the 2048 waveform samples into an appropriate cell. The results, 
divided by the total number of samples, is the amplitude probability function, examples 
of which are given here. The computer program determines the mean value of the fluc- 
tuation amplitudes and shifts the zero of amplitude to the nearest amplitude cell for 
plotting. The p'' moment of f(x) is defined as 


ct(p) s 


/■ 


x%(x)dx 



lp/2 
x^f (x)dx 


( 1 ) 
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where the mean value of the variable x is zero. The computer program calculates mo- 
ments up to p = 7. The moments for some common probability density fvinctions and 
for some common waveforms are given in table n. The standard deviation a(2) is a 
measure of the width of the distribution function, or the degree of spread about the mean 
value. The skewness a(3) is a measure of the asymmetry of the distribution function. 
Negative values of skewness imply distribution functions that are lopsided to the right 
of the mean value; positive values imply the opposite. The kurtosis ct( 4) is a measure 
of how much time the signal spends in the wings of the distribution: A kurtosis of 1. 0 


TABLE n. - MOMENTS FOR SOME COMMON PROBABILITY 


DENSITY FUNCTIONS AND WAVEFORMS 


Moment, 

P 

Square wave 
of peak am- 
plitude A 

Sinusoid of 
peak ampli- 
tude A 

Full -wave 
rectified sine 
wave of peak 
amplitude A 

Gaussian 

random 

noise 

Moment 

1 - Mean 

0 

0 

0. 637 A 

0 

2 - Standard deviation 

A 

0. 707 A 

0. 308 A 

a 

3 - Skewness 

0 

0 

-0. 497 

0 

4 - Kurtosis 

1.0 

1.5 

1.93 

3 

5 - Superskewness 

0 

0 

-2. 11 

0 

6 - Superkurtosis 

1.0 

2.5 

5.07 

15 

7 - Hyperskewness* 

0 

0 

-7.69 

0 

8 - Hyperkurtosis 

1.0 

4. 375 

15.8 

105 


implies that the signal spends all its time at the extreme values of the range (in the 
manner of a square wave); higher pos^itive values imply increasingly centralized, 
delta-function-like distributions. The higher moments, above p = 4, do not appear to 
add new physical content to the lower moments, and therefore will not be considered 
further herein. The moments for a square wave, for a sinusoid, for a full-wave recti- 
fied sine wave, and for Gaussian random noise are shown as examples in table II. 

The moments of the probability distributions for the runs listed in table I are 
given in table m. The standard deviation is given directly in volts and is (the first 
entry in table m) as large as 45 percent of the electrode voltage. Most deviations 
range from 10 to 30 percent of the electrode voltage. These very violent potential 
fluctuations must be expected to play a role in the stability and transport of the plasma. 
The kurtosis of nearly all runs is close to 3. 0, which is characteristic of Gaussian 
random noise, for both potential and density fluctuations. An exception is the first run, 
in which the kurtosis is close to that of a sinusoid. As shown later, this run is char- 
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TABLE m. - MOMENTS OF PROBABILITY DISTRIBUTIONS 


Run 

Figure 

Capacitive-probe potential fluctuations 

Langmuir^probe density 






- 



fluctuations (channel 3) 



Near Langmuir probe 

30^ above Langmuir probe 






(channel 1) 

(channel 2) 

Standard 

Skew- 

Kurtosis 









deviation, 

ness 




Standard 

Skew- 

Kurtosis 

Standard 

Skew- 

Kurtosis 

V 





deviation, 

ness 


deviation, 

ness 







V 



V 






1 

— 

780 

-0. 12 

1.6 

900 

-0.08 

1.6 

3.2 

-0. 16 

2.3 

2 

11 

270 

-.47 

2.2 

300 

-.49 

2. 2 

1.2 

-.03 

3.7 

3 

12 

no 

-.53 

3.4 

160 

.02 

2.8 

2.7 

. 21 

3.2 

4 

— 

100 

-.52 

3.2 

no 

-. 53 

3.3 

.21 

. 22 

3, 1 

5 

— 

210 

.27 

2.8 

230 

. 55 

3.7 

4. 1 

. 55 

3.9 

6 

— 

62 

-.04 

3.0 

no 

. 11 

2.6 

8.5 

.72 

3. 1 

7 

— 

99 

11 

2.7 

120 

-.01 

2.7 

16 

.97 

4.0 

8 

14 

95 

-. 18 

3.0 

81 

-.43 

3.5 

2. 1 

.24 

3.2 

9 

15 

no 

-.25 

2.4 

140 

-. 21 

2.5 

5.2 

,94 

4.4 


acterized by an extremely coherent rotatir^ spoke at a single frequency. For the ap- 
proximately 30 experimental runs analyzed, the skewness of the density fluctuation 
distribution (channel 3 in table III) seems to assxmie high positive values, but that of the 
potential fluctuations tends to be smaller and to depend on the circumstances of the run. 


Low-Pressure Mode 

Some examples of the distribution functions listed in table HI are given for the low- 
pressure mode of operation. Figure 11 gives the distribution functions for the second 
run in tables I and III. Each point plotted represents the fraction of the 2048 sample 
points in the cell that corresponds to the potential or density amplitude plotted on the 
abscissa. This run illustrates a condition in which the potential distribution is mod- 
erately highly skewed, as a result of the raw potential waveform being cycloidal with 
the cusps pointing downward. The density distribution function is nearly symmetric. 
Figure 12 gives the distribution fimctions for the third run in tables I and IE. This run 
illustrates a condition for which the kurtosis of all the distributions is close to that of 
Gaussian random noise but for which the skewness of the two potential distributions 
differs substantially. Finally, figure 13 gives the distribution functions for the ninth 
run in tables I and III. This run illustrates a distribution for which the density fluc- 
tuations are highly skewed but the potential fluctuations much less so. 
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High-Pressure Mode 


A characteristic example of the amplitude probability distribution function for the 
high-pressure mode (HPM) of operation (eighth entry in table HI) is given in figure 14. 
The distribution functions of the density and potential fluctuations tend to have a 
kurtosis near 3. 0, which is characteristic of Gaussian random noise, but they also 
have nonzero values of skewness with the potential distributions beii^ negative and 




Figure 14. - Probability density function for run 8 of table I. 
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skewed to the right and density distributions beii^ positive and skewed to the left. The 
tendency of both modes of operation to have nonzero values of skewness can be under- 
stood as a result of the cycloidal, cusplike waveforms associated with the motion of a 
rotating spoke as it moves past the probes (refs. 2 and 3). The tendency of the skew- 
ness to have opposite signs for the density and potential fluctuations can be understood 
if both waveforms are cusplike but point up or down in opposite directions. 


AZIMUTHAL DISPERSIVE CHARACTERBTICS OF PLASMA 
Theoretical Considerations 

Low-frequency potential fluctuation spectra were obtained by processing the signals 
^j^(t) and from the two capacitive probes in channel 1 and channel 2, respectively. 

A schematic of the data-processii^ system is shown in figure 9. The computer output 
consists of plots of the autopower spectnun, the cross-power spectrum (both amplitude 
and phase), and the squared coherency spectrum between channels 1 and 2. 

The underlying principles behind fast- Fourier -transform analysis of plasma fluc- 
tuation data are described in detail in references 5 and 21 to 23. Here we give a brief 
review of digital spectral analysis and how it may be used to analyze and interpret 
plasma fluctuation data. The approach is to take the two capacitive probe signals 
^j^(t) and digitize them, and to compute their respective Fourier transforms 

$i(co) and by using the fast- Fourier-transf orm algorithm. The autopower spec- 

trum^ of channel 1 is given by Pj^^(co) = $ j(uj)$ j(o)), where the asterisk denotes a com- 
plex conjugate. The autopower spectrum for channel 2 is computed in a similar way. 

Of particular importance is the cross-power spectrum between ^j(t) and ^2^ 
since it can suggest the presence of propagating waves. The cross-power spectrum is 
given by ^ complex quantity. In this transport work, de- 

scribed later in this section, we considered the real and imaginary parts of the cross- 
power spectrum, which are known as the co- and quadrature spectra, respectively. 
However, in studying the properties of the fluctuations, it is convenient to express the 
cross-power spectrum in terms of an amplitude spectrum 1P|2(^)| ^ phase 

spectrum. From the definition of = ^i(to)^^(co), it is clear that the cross- 
amplitude spectrum is given by | Pi2(“^) 1 = 1 I phase spectrum by 

012(co) = 0j^(w) - where 0j(u>) is the phase of ^j(co) and 02^*^) phase of 

^Strictly speaking, the autopower spectrum is defined as 
Pj^j(w) = limZX E j^(w)<5^(aj)]|, where T is the time duration of the ^j^(t) waveform 

used to compute <^1(0)). A similar definition applies to the cross-power spectra. 
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# 2 (^)* We use the important fact that the phase spectrum preserves phase information 
in the form of a phase difference to determine the wave number of each wave in the 
measured fluctuation spectrum. For example, if a wave is present at frequency co, 
the phase shift that the wave undergoes in propagating from one probe to another will 
be given by 0^2^(jo) and will be equal to k(co)As; thus = k(co)As. This equation 

suggests then that k{o)) can be determined by takii^ the computer-generated phase 
spectrum 0 dividing by As, the known distance of separation between the 

probes. In cylindrical geometry, and for probes displaced azimuthally by A0, 

As = Rp A9 and = m(o))/Rp, where R^ is the radial position of the probes and 
m(co) is the azimuthal mode number of the wave. With this information it is easy to 
show that m{o)) = 0 j2(^^)/A0. 

To summarize, from the peaks in either the autopower spectrum or the cross- 
amplitude spectrum amplitude and frequency can be determined for each 

of several waves in the fluctuation spectrum. From the phase spectrum 0 ^ke 
corresponding wave number k(a;) = or azimuthal mode number 

m(co) = 0^2^.^)/^^ can be determined. Knowledge of the wave number permits com- 
putation of the phase velocity co/k at each frequency. The phase velocity can be 
written 


and can be obtained from the computer-generated phase plots of versus co. 

From the magnitude, direction, and parametric variation of the phase velocity the 
wave can often be identified. 

Another spectral quantity of some importance is the coherence spectrum, defined 
as “ I Its properties are discussed in both ref- 

erences 23 and 24. Basically, it is a measure of the degree of mutual coherence be- 
tween the signals (pj^(t) and cp^t)^ with the difference in phases ignored. The coher- 
ence spectrum is bounded by zero (incoherent case) and unity (perfectly coherent 
fluctuations). Values between these two ideal limits correspond to a partially coherent 
case. In the following discussion, computer -generated plots of | shown. 

High values of coherence are indicative of high signal-noise ratios, and conversely. 


Low-Pressure Mode 


Some characteristics of the peaks in the autopower spectra of the experimental 
runs listed in table I are shown in table IV. Some had more than one peak in the spec- 
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TABLE IV. - CHARACTERISTICS OF TRANSPORT SPECTRAL PEAKS 


Run 

Figure 

Peak 

fre- 

quency, 

kHz 

Azimuthal 
phase ve- 
locity at 
peak, 

cm/sec 

Channel 1 
autopower 
peak am- 
plitude, 

V 

Channel 2 
autopower 
peak am- 
plitude, 

V 

Channel 3 
autopower 
peak am- 
plitude, 
m“^ 

Peak transport 
rate, 

number/m^‘ sec 

Total cumulative 
transport rate, 
number/m^ • sec 

Total cur- 
rent (cal- 
culated), 
A 

Direction 
of radial 
transport 

1 

21 

20 

1. 6x10® 

350 

400 

1. 2x10^^ 

37x10^® 

15X10^'^ 

0.53 

Inward 

2 

20 

15 

.61 

190 

210 

.56 

22 

4.2 

. 15 



3 

25 

42 

«0 

24 

22 

.49 

«0 

.08 

.003 



4 

29 

— 

«0 

--- 

... 



.027 

.001 



5 

--- 

30 

2.9 

57 

65 

.60 

1.6 

1 ] .68 

.024 





180 

1.4 

18 

20 

.53 

1.6 

J 




6 

27 

130 

2. 1 

12 

20 

1.4 

12 

\ . 26 

. 93 





180 

2. 1 

18 

42 

2.4 

48 

J 




7 

— - 

140 

4.0 

29 

48 

5.8 

47 

\ .33 

1. 18 





194 

2.0 

22 

24 

3.8 

72 




f 

8 

26 

145 

1.3 

13 

8.8 

.31 

-3.4 

1-2.7 

.097 

Outward 

1 


195 

1.7 

7 

13 

.51 

-2.5 

J 




9 

23 

11 

. 56 

54 

69 

.63 

-3.5 







121 

-3.4 

9 

21 

1.2 

1. 1 

^-3.7 

. 13 

Outward 



175 

-1.5 

4.9 

7.3 

.90 

-4.6 

J 





trum. The peak frequency and azimuthal phase velocity (calculated from eq. (2)) are 
listed for each peak in columns 3 and 4. The next two columns list the root-mean- 
square (rms) amplitude of each peak in volts. The rms amplitude of the peaks in the 
density spectrum (column 7) are given in units of 10^*^ per cubic meter. The amplitudes 
of the potential peaks range from about 10 volts to several hundred volts. Run 4 is an 
example of a spectrum so turbulent that no separate peaks could be identified. The 
specific figures illustrating the entries in table TV are identified in the table. 

The fluctuation spectrum in the plasma with electrodes 5 and 6 in place, an anode 
voltage of 2 kilovolts, a neutral background density of 3. 5x10^^ particles per cubic 
centimeter (1. 06x10"^ torr) - run 5 in tables I and IV - is shown in figure 15. The peaks 
in the cross-power spectrum that have high coherency are at 30 and 180 kilohertz. The 
mode at 30 kilohertz propagates with an azimuthal phase velocity of 3. 9x10” centime- 
ters per second in the direction of the E/B drift with a radially outward-pointing elec- 
trie field. The wave at 180 kilohertz propagates with an azimuthal velocity of 1. 4x10 
centimeters per second in the same direction. The spectral band from 0 to 75 kilohertz 
has high coherency and very little jitter in the phase plot, implyii^ the presence of a 
continuum of eigenmodes. The continuum propagates along with the 30-kilohertz mode. 

It has an amplitude of about 15 volts rms. The dominant mode at 30 kilohertz is oscil- 
lating with an amplitude of 60 volts rms. The power spectriun falls off gradually at 
250 kilohertz, where it corresponds to a noise level of about 4. 5 volts rms. The radial 
profile of the floating potential measured with a Langmuir probe indicates that the plasma 
boundary is subjected to a radially outward-pointii^ electric field of no less than 90 volts 
per centimeter. For a magnetic field of 0. 67 tesla at the probe, this corresponds to an 
E/B drift velocity of 1. 3x10 centimeters per second. The magnetic field in the bumpy 
torus is highly nonuniform both radially and axially. The region in which the probe was 
located had the lowest magnetic field in the entire containment voliune. Any higher mag- 
netic field in the plasma volume will give lower values of the E/B drift velocity, given 
the electric field measured at the probe location. In all calculations, a magnetic field 
of 0. 67 tesla, corresponding to the region where the floating potential profile and the 
wave dispersion were measured, was used. 

The wave dispersion measured when the sector diametrically opposite to the one 
containii^ the probe assembly had an electrode ring (electrode in section 10 only) is 
shown in figure 16 and corresponds to run 7 in tables I and IV. The coherent modes 
occur at 140 and 194 kilohertz and propagate with velocities of 4. OxlO” and 2. OxlO” 
centimeters per second, respectively, in the direction of the E/B drift velocity. The 
radial electric field for this case was 110 volts per centimeter and points radially out- 
ward. It corresponds to an E/B drift velocity of 1. 6x10® centimeters per second. 

Figure 17 shows the power spectra of the fluctuations present in the plasma gener- 
ated with electrodes 1 and 7 in place at a background neutral density of 7x10^^ particles 
per cubic centimeter (2. 1x10“® torr) and an anode voltage of 1 kilovolt (run 3 in tables 
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Frequency, kHz 

Figure 16. - Azimuthal wave dispersion for run 7 of table L 
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Figure 17. - Azimuthal wave dispersion for run 3 of 
table I. 


I and IV). The strongest peak with an amplitude of 25 volts rms occurs at 42 kilohertz. 
The phase spectrum indicates that k^ = 0 at 42 kilohertz. The two other coherent oscil- 
lations present are at 72 and 84 kilohertz. The mode at 72 kilohertz propages with a 
velocity of 5. 5x10^ centimeters per second in the direction of the E/B drift with a ra- 
dially outward -pointing electric field. Low coherency in the spectral band 0 to 35 kilo- 
hertz indicates the presence of noise. A measured radial electric field of 120 volts 
per centimeter corresponds to an E/B drift velocity of 1. 8x10° centimeters per second. 
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High-Pressure Mode 


The wave dispersion when all the electrodes except sector 4 were used to generate 
the plasma is shown in figure 18 (rim 4 in tables I and IV). The coherency plot shows 
that the spectral band 1 to 10 kilohertz has high coherency, and the phase spectrum 
shows that k^ = 0 in this band. The cross-power spectrum lacks any peaks, though 
the amplitude of the fluctuations in the 1- to 10-kilohertz band is about 15 volts rms. 




Figure 18. - Azimuthal wave dispersion for run 4 of 
table I. 
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The power spectrxim falls off gradually between 10 and 50 kilohertz. At 50 kilohertz, 
the amplitude is about 4 volts rms, the coherency is low, and the phase is jittery - 
characteristics of noise. Examination of a wider spectral band shows that figure 18 is 
representative of the fluctuation spectriun with no additional peaks at higher frequencies 
than those shown. A dc radial electric field of 80 volts per centimeter was present at 
the probe location, yielding an E/B drift velocity of 1. 2x10 centimeters per second. 

Figure 19 shows a case like that of figure 17 except that the neutral background 
density has been changed to 1. 6X10'*’ particles per cubic centimeter (4. 8x10"^ torr), 
run 8 in tables I and IV. The autopower spectrum of channel 1 has peaks at 50, 145, 
and 195 kilohertz. The mode at 50 kilohertz is characterized by k„ = 0. The waves 
at 145 and 195 kilohertz are propagating with phase velocities of 1. 3xl0” and 1. 7x10° 
centimeters per second in the direction of the E/B drift with a radially outward-pointing 
electric field. 


FLUCTUATION -INDUCED PARTICLE TRANSPORT 
Theoretical Considerations 

The key ideas underlying our approach to measuring fluctuation-induced transport 
are described in detail in reference 25. Low-frequency (w « fluctuations are as- 
sumed so that a particle’s fluctuating velocity may be modeled in terms of an E,/B drift, 
where Ej^ is a fluctuating electric field and B is the static, confining, toroidal mag- 
netic field. The time-averaged particle flux is then given by 

(nE|) 1 


where C_t;,( 0) is the cross-correlation function between density and electric field flue- 
tuations evaluated for zero time delay. Our goal is to relate {nv; to experimentally ob- 
servable spectral properties of the fluctuations. To do this we recall that the cross- 
correlation function between density and potential fluctuations may be expressed in terms 
of the inverse Fourier transform of the cross-power spectrum between density and elec- 
tric field fluctuations Pj^gCco): 


/ OO 




(4) 
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Substituting equation (4) into equation (3) yields 

<■;?) i: PnE(w)dw 


( 5 ) 


where the time-averaged particle flux is expressed as an integral over all frequencies 
of the real part of the cross-power spectrum between density and electric field fluctua- 
tions. 

In the NASA Lewis bumpy-torus experiment we have the capability to measure po- 
tential rather than electric field fluctuations. To the extent that the quasi-static approxi- 
mation holds, the cross-power spectrum between density and electrical 

field fluctuations “lay be expressed in terms of the cross-power spectrum be- 
tween density and potential fluctuations follows: 

where kj^(co) is the perpendicular wave number. 

Substituting equation (6) back into equation (5) yields the final result 

<nv) T(w)dw (7a) 

which states that the time-averaged particle flux is given by an integral over all fre- 
quencies of the transport spectral density function T(w). The quantity T(w) is defined 
as 


T(co)= -?kj^(co)Q (co) (7b) 

B ^ 

where the quadrature spectrum Qjj^(w) is defined as the imaginary part of Pji^(<^)> The 
factor 2 in equation (7b) takes into account the "negative frequency" contributions to the 
particle transport. (T(w) is an even function of frequency. ) The function T(o>) is a 
"transport spectral density function" in that its dimensions are number of particles per 
unit time, per unit area, and per unit radian frequency. Thus a visual inspection of a 
plot of T(w) not only will indicate which parts of a fluctuation spectrum are contributing 
to transport, but also will indicate in a quantitative way how much transport is associated 
with each spectral band. 

Our approach in generatii^ the transport spectral density function is as follows: 

First, the signals (channel 1 and channel 2) from the two capacitive probes displaced 30° 
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apart in azimuth are digitized and the perpendicular wave number kj^(oj) is determined 
by using the procedures reviewed in the section AZIMUTHAL DISPERSIVE CHARAC- 
TERISTICS OF PLASMA, Next the signal from the Langmuir probe (channel 3) is digi- 
tized and its Fourier transform computed. From the Fourier transforms of the density 
fluctuation and the potential fluctuation signal in channel 1, the quadrature spectrum 
Q (w) is computed by using the procedures described in reference 23. The transport 
spectrum is then generated by carrying out the multiplication indicated in equation (7b). 

Although equation (7b) is useful in terms of how to compute a transport spectrum, 
an alternative formulation is advantageous in terms of understanding which effects result 
in large or small amounts of transport. For example, it is shown in reference 25 that 
the transport due to a small band of frequencies centered at w is given by 

k,(a)) 

T(o)) = — - I 

Therefore, the transport associated with this band of frequencies centered about w de- 
pends on the following factors: the product of the rms values of density and potential 
fluctuations, the sine of the phase angle Q^n^Cw) between the density and potential fluc- 
tuations, and the degree of mutual coherence I between the density and potential 

fluctuations in the spectral band under consideration. The factor k (oi) appears since 

I I 

I — rms^^^^’ model assumes v = E^B. 

In subsequent figures, computer-generated plots of transport spectra are shown in 

O 

the lower right corner. Either autopower spectrum of channel 3, or 

n(<x))^(u)) |yjj^(o)) I, the cross-power spectrum between channels 1 and 3, is plotted in 
the upper left corner. The phase-difference spectrum 0 !j^^(w) for density and potential 
fluctuations is plotted in the lower left corner and corresponds to the phase spectrum 
between channels 1 and 3. Since the transport depends on sin Oj^^(o)), values of o'j^^(o>) 
near ±90° are optimum in terms of transport. Finally, the square of the degree of mutual 
coherence between density and potential fluctuations \ is plotted in the upper 

right corner and corresponds to the square of the degree of mutual coherence between 
channels 1 and 3. As explained in reference 25, the degree of mutual coherence between 
density and potential fluctuations enters into equation (8) because we are considering a 
polychromatic case rather than a monochromatic case. In the monochromatic case the 


n 

Since the data have been digitized, we actually plot a transport spectrum 
(Number of particles/ Area- Time) not a transport spectral density function. To deter- 
mine the total transport associated with a given fluctuation spectrum, we sum, rather 
than integrate, the transport spectrum over all frequency values. 
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3 

degree of mutual coherence would be unity. 

The transport spectrum T(o)) is a real quantity and may be either positive or nega- 
tive depending on whether the transport is in the inward or outward direction. Equa- 
tion (8) indicates that only two quantities determine the s^n of T(co) and hence the direc- 
tion of transport: The first quantity is the phase difference between density and 

potential fluctuations. If changes sign, T(w) changes sign since T(a>) is propor- 
tional to sin . The second quantity determinii^ the sign of T(o)) is the sign of 

kj^(o)), which in effect determines the direction of the fluctuating electric field Ej^ in 
equation (3). 


Low-Pressure Mode 

The transport spectrum when all the electrodes except that in sector 4 were in 
place, and with the electrodes biased to 1 kilovolt at a neutral background density of 
7x10^^ particles per cubic centimeter (2. 1x10”^ torr), is shown in figure 20 (run 2 in 
tables I and IV). This figure is a plot of the cross-power spectrum, squared coherency 
spectrvun, phase spectrum, and transport spectrum obtained from the probes in chan- 
nels 1 and 3. The transport is radially inward in this case. The cross-power spec- 
trum is dominated by a peak at 15 kilohertz that corresponds well with a similar peak 
in the coherency spectrvun. This implies that there is a coherent wave mode in the 
potential and density fluctuations. The phase spectrum reveals that at 15 kilohertz the 
phase difference between density and potential fluctuations is 120®. As the transport 
spectrvun indicates, there is negligible transport in the whole spectral band outside 
the peak at 15 kilohertz. Keeping all other parameters the same, we changed the anode 
voltage to 2 kilovolts (rvui 1 in tables I and IV). Figure 21 shows the corresponding 
transport spectrvun, which is also radially inward. The transport spectrvun and the 
cross -power spectrvun show that the coherent wave that causes large amounts of trans- 
port now occurs at 20 kilohertz. The phase difference between density and potential 
fluctuations is now 40® and the amplitude of the peak in the cross-power spectrvun is 
again 22 decibels above the noise level. The spectral band from 30 to 100 kilohertz 
has low coherency and a jittery phase spectrvun - characteristics of noise. In this 
frequency interval, the transport spectrum indicates very little transport of plasma. 
Figure 22 shows the dispersion spectra for the conditions corresponding to figure 21. 


^Other properties of coherency spectra are reviewed in ref. 24. In particular, low 
values of |yjj^(w)| i»ay be due to noise in the system, a departure from a linear rela- 
tion between density and potential fluctuations, and/or so-called multiple inputs. 
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Phase spectrum, rad Cross-power spectrum, n(o;)$(cj)|yp^(£»;)|. 



Frequency, kHz 

Figure 20. - Transport spectrum for run 2oftableI. 
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These indicate a very coherent wave propagating at 1. 6x10 centimeters per second 
at 20 kilohertz. 

To study the transport with different electrode ring configurations, we kept the 
probe assembly fixed in sector 4 and changed the location of the electrodes from sec- 
tors 3 and 9 (figs. 23 and 24 and run 9 in tables I and IV) to sectors 1 and 7 (fig. 25 
and run 3 in tables I and IV). The neutral background density was kept at 7x10^ ^ par- 
ticles per cubic centimeter (2. 1x10”^ torr) and the anode voltage, at 1 kilovolt. 'When 
the nearest electrode is one sector away, the cross spectrum (fig. 23) has peaks at 
121 and 178 kilohertz, which correspond to high coherency in the squared coherency 
spectrum. There is hardly any transport at 121 kilohertz because of the small phase 
difference between the potential and density fluctuations. Figure 24 shows the dispei — 
Sion spectra for the conditions in figure 23. These spectra indicate a relatively co — 

f* 

herent disturbance propagating at 0. 56x10 centimeters per second at 11 kilohertz. 

The net transport was radially inward in this case. The transport spectrum when the 
electrodes are three sectors away on each side of the probe is shown in figure 25. The 
net transport at the probe location was outward in this case. The autopower spectrum 






Rgure 23. - Transport spectrum for run 9 of table I. 
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of the ion density fluctuations has a dominant peak only at 42 kilohertz. However, the 
capacitive probe at the same point in the plasma picks up signals at 42, 72, and 84 
kilohertz, as shown in figure 17. The disturbance at 42 kilohertz is characterized by 
kg = 0, since the value of approximately zero at this frequency (fig. 17). It 

is immediately noticeable that an approximately zero kg in the transport equation 
(eq. (7b)) results in very low transport at 42 kilohertz in figure 25. 


High-Pressure Mode 

The transport spectrum in the high-pressure mode for the conditions of figure 25 
is shown in figure 26 (run 8 in tables I and IV). The density fluctuation spectrum has 
peaks at 145 and 195 kilohertz. The net radial transport is outward in this case. Even 
though in figure 26 the amplitude of the disturbance at 145 kilohertz is smaller than that 
at 195 kilohertz, the absolute amount of transport is more at 145 kilohertz. The reason 
for this is that at 145 kilohertz the phase difference between density and potential is 
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Frequency. kHz 

Figure 26. - Transport spectrum for run 8 of table I. 

optimum (90°) for maximum transport. The transport associated with the 145- and 
195-kilohertz peaks is negative since, at these frequencies, the phase difference be- 
tween density and potential fluctuations is negative. 

Figure 27 shows the transport spectrum present in the plasma when the plasma 
was generated by an electrode in a sector (10) diametrically opposite the probe as- 
sembly. Disturbances that have peaks in the cross-power spectrum and high coher- 
ency are at 130 and 180 kilohertz. The dominant transport occurs at 180 kilohertz, and 
the net transport is radially inward. The phase difference between density fluctuations 
and potential fluctuations is an optimum of 90°. Figure 28 shows the dispersion spectra 
for the conditions of figure 27. These spectra indicate that virtually all disturbances 
up to 250 kilohertz are propagating at 2. 1x10 centimeters per second. 

Figure 29 illustrates the fourth run in tables I and IV, the dispersive characteris- 
tics of which are shown in figure 18. In this case, the coherency is low and the trans- 
port is very low in magnitude and quite turbulent in the sense of having no dominant 
peaks in the spectrum. A summation over all frequency bands of the transport spec- 
trmn yields a net inward transport of particles. 
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Total Ion Transport Rates 

The last four columns of table IV contain information about the absolute rates of 
transport as measured with this technique. The first of these (colmnn 8) contains the 
transport rate associated with each of the dominant peaks that are evident in the trans- 
port spectrum. In the third and fourth runs, no peaks appear above the background tur- 
bulent noise, and in the last five runs there is more than one peak in the spectrum. 
Positive values of transport are radially inward, and negative values are radially out- 
ward. The total cumulative transport summed over all frequency bands of the spectrum 
(column 9) represents the flux of ions at the probe. In collecting the raw data, care was 
taken to assure that there was no significant fluctuation activity at higher frequencies 
that might have contributed to the transport. The last column gives the direction of the 
total transport at the probe. 

Charged particles are lost from the plasma at a total rate equal to the current drawn 
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by the power supply since previous work has shown that no ambipolar currents flow to 
the electrodes (ref. 4). This total current can be compared with the total, cumulative 
transport rate only if it is assumed that the probes were located at a typical position. 
This assumption is not valid since the magnetic field at the probe was weaker than any- 
where else in the plasma volume, the probes were inside the plasma boundary rather 
than on or outside it, and the magnetic field and plasma properties vary both around the 
minor circumference of the plasma and axially along the magnetic field lines. 

Although these factors mitigate against a comparison of the measured transport 
rates with the electrode current, it is nonetheless Interesting to do so. The surface 
area of the bumpy-torus plasma is 2. 23 square meters. If this area is used to convert 
the trainsport rates in column 9 of table IV to current flows, the currents shown in 
column 10 result. Figure 30 plots the measured electrode current from table I as a 
function of the calculated transport current from table IV. The transport currents are 
generally within an order of magnitude of the measured currents, about as good agree- 
ment as can be expected in view of the nonuniformity of the plasma about the major and 
minor circumferences of the plasma. 



Figure 30. - Measured transport current as a function of electrode current flowing 
to power supply for experimental runs tabulated in table I. 
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DISCUSSION 


This report has described the first application of a technique by which the radial 
particle flux in a hot-ion, steady-state plasma can be measured directly. These data 
were exploratory and no attempt was made to obtain parametric variations of the dis- 
persion or transport characteristics of the plasma. A wide range of phenomena have 
been revealed, including rotating disturbances that move with velocities comparable to 
the E/B drift velocity, transport that is radially inward or outward, and radial transport 
that is associated with random turbulence as well as with strong coherent peaks in the 
fluctuation spectra. 

The probability density fxmctions of the fluctuating density and potential amplitude 
are valuable as crosschecks for signal saturation of the raw data. When s^nal satura- 
tion occurs, it is very evident from the large, delta-function-like wings at the extreme 
edges of the probability density function. The distribution functions of both density and 
potential generally had a kurtosis of about 3. 0, which is characteristic of Gaussian ran- 
dom noise, and nonzero values of skewness, which result from the cusplike waveforms 
that have been observed in this plasma (refs. 2 and 3). (Random Gausian waveforms 
have zero skewness. ) The first run in table HI is an example in which the raw wave- 
form was very coherent and periodic, with essentially zero skewness (a sinusoidal 
signal) and kurtosis appropriate to a sinusoid. 

When the operatii^ conditions were varied, the plasma exhibited electrostatic po- 
tential fluctuations (1) at a set of discrete frequencies, (2) at a continmun of frequencies, 
and (3) as incoherent turbulence. The coherent discrete modes propagated azimuthally 
in the direction of the E/B drift, with an electric field that usually pointed radially out- 
ward between the plasma and the magnet dewars. Electric fields in the core of the 
plasma could not be measured with Langmuir probes. A nonperturbing technique, such 
as ion beam probing, may be required to yield the potential profile in the core of the 
plasma. 

In several runs (of which figs. 15 and 19 are characteristic examples), coherent, 
discrete modes with k^ = 0 have been observed. In such cases, there is a very large 
peak in the low-frequency potential fluctuation spectrum (below 100 kHz), but there is 
essentially no time delay between the fluctuations at the site of the two probes. This 
oscillation of the plasma in unison may be a manifestation of the continuity- equation 
oscillation, which has been observed in a similar plasma (ref. 26). The possibility that 
the probes are seeing a radially or axially propagating wave cannot be ruled out on the 
basis of the data taken for this report, but previous work (ref. 3) demonstrates that co- 
herent low-frequency phenomena, in phase around the major circumference of the 
torus, do occur. A measurement of k|| and k^, will be helpful in confirming the pres- 
ence of the continuity -equation oscillation in this plasma. 
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The data on the particle transport spectrum revealed a wide variety of phenomena: 

(1) All transport occurring in a single, sharp peak representing radially inward 

transport (figs. 20 and 21) 

(2) Transport occurring in two or more peaks, one carrying particles inward and 

one outward (fig. 23) 

(3) Spectra with two peaks, both carrying particles inward (fig. 27) 

(4) Spectra with two peaks both carrying particles outward (fig. 26) 

(5) Spectra with no prominent peaks, in which the total transport is small and caused 

by a continuum of turbulent fluctuations (figs. 25 and 29) 

Particle transport was radially inward for a rnajority of the cases examined, an 
interesting result. This inward transport could be caused by the stroi^ radial and 
azimuthal electric fields associated with the rotatii^ spokes in the plasma (ref. 3) and 
may be prevalent in the sector in which the probe was located because there was no 
electrode rit^ in that sector, for the data discussed in this report. The transport calcu- 
lated from the model of reference 25 appears to account for substantially all the particle 
flux intercepted by the anode rings and the walls. The data in figure 30 indicate that the 
magnitude of the transport currents, if the probe position is typical of the entire surface 
of the plasma, is within order-of -magnitude agreement with the currents actually flow- 
irg to the power supply. 


CONCLUSIONS 

Simultaneous density and potential fluctuation data have been taken in a hot-ion, 
steady-state, bumpy-torus plasma subject to strong, crossed electric and magnetic 
fields. The statistical properties of the density and potential fluctuations, the disper- 
sion of the propagating disturbances in the electrostatic potential, and the radial trans- 
port spectrum have been calculated. The transport spectrum is based on the model pre- 
sented by Powers, according to which all radial transport is due to E/B drift caused by 
potential fluctuations and the resulting electric fields. All measured and calculated 
quantities, includii^ the potential and density amplitudes and the transport rate, were 
obtained in absolute units. ' 

At low frequencies (below 100 kHz), a nonpropagating azimuthal mode was observed 
with kg = 0. This may be related to the continuity- equation oscillation, which is known 
to have kg = 0. Additional measurements are needed of k^. and kjj to unambiguously 
identify this azimuthal mode. The azimuthally rotating spokes had velocities consistent 
with the E/B drift velocity at the probe. These rotating spokes were often associated 
with the peaks in the radial particle transport spectrum. The peaks for which kg = 0 
did not have significant transport associated with them. 
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The most surprising feature of the radial transport of ions is that much of it was 
radially inward in this plasma. The model of Powers appears to account for all the 
transport in this plasma, because the calculated radial transport rates are in order- 
of -magnitude agreement with the electrode currents actually drawn from the plasma. 

The highest rates of radial transport were associated with sharp, coherent peaks in the 
transport spectrum, and the transport associated with random, turbulent spectra with- 
out peaks was much smaller in magnitude. 

The exploratory data reported herein have proven the basic value of the radial 
transport model given by equation (3) and of the hardware and software systems devel- 
oped to acquire and analyze the required data. Future work along these lines should in- 
clude a study of parametric variations, includit^ the transport rate as a function of elec- 
trode voltage, plasma density, and other independent variables that affect the plasma 
characteristics. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 4, 1978, 

506-25. 
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APPENDIX - SYMBOLS 


B 

C 

E 

E(A, B) 
e 
f 
k 

m(co) 

n 

p, ,(<..) 

Plj(cc) 

Q(w) 

R 

AR 

T 

T(w) 

t 

V 

V 
X 


a 


Ad 


\ 

<y(p) 


T 

$ 


magnetic field strength, T 
cross- correlation fimction, eq. (3) 
electric field strength, V/m 
expected value of A, B 
electronic charge, C 
amplitude probability function 
wave number, m”^ 
mode spectrum of propagating waves 

_3 

number density, cm 
autopower spectrum 
cross-power spectrum 
quadrature spectrum 
radius from plasma axis, m 
distance between probes, m 
kinetic temperature, eV 
transport spectral density function 
time, sec 
voltage, V 
velocity, m/sec 
fluctuation amplitude 
phase angle, rad 

coherence spectrum of signals i, j 
phase spectrum 

angular distance between probes 
wavelength, m 

**p*’ " moment of distribution fxmction, eq. (1) 

time interval, sec 

Fourier transform of potential, V 


<p potential, V 

CO frequency, rad/ sec 

Subscripts: 
a anode 

ci ion cyclotron 

e electrons 

i ions 

n particle number density 

p probe 

p0 phase velocity in azimuthal direction 

r radial direction 

rms root mean square 

9 azimuthal direction 

(p potential 

II parallel to magnetic field 

i perpendicular to magnetic field 

1 capacitive probe, channel 1 

2 capacitive probe, channel 2 

3 Langmuir probe, channel 3 
Superscripts: 

“ time-averaged value 

~ fluctuating value 

* complex conjugate 
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